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Abstract 

We discuss orbital entanglement in mesoscopic conductors, focusing on the effect of dephasing. The 
entanglement is detected via violation of a Bell Inequality formulated in terms of zero-frequency current 
correlations. Following closely the recent work by Samuelsson, Sukhorukov and Biittiker[l], we investigate 
how the dephasing affects the possibility to violate the Bell Inequality and how system parameters can 
be adjusted for optimal violation. 
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1. Introduction 

Entanglement is one of the most intriguing features predicted by quantum theory [2, 3, 4]. It leads to 
correlations between distant particles, which can not be described by any local, realistic theory [5]. This 
nonlocal property of entanglement has been demonstrated convincingly in optics [6, 7], where entangled 
pairs of photons have been studied over several decades [8]. Apart from the fundamental aspects, there is a 
growing interest in using the properties of entangled particles for quantum cryptography [9] and quantum 
computation [10]. 

Recently, much interest has been shown for entanglement of electrons in solid state systems. A controlled 
generation and manipulation of electronic entanglement is of importance for a large scale implementation 
of quantum information and computation schemes. Electrons are however, in contrast to photons, massive 
and electrically charged particles, which raises new fundamental questions and new experimental challenges. 
Most existing suggestions [11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21] are based on creating, manipulating and 
detecting spin-entangled pairs of electrons. This requires experimental control of individual spins via spin 
filters or locally directed magnetic fields on a mesoscopic scale. 

To overcome these difficulties, very recently spin-independent schemes for creating and detecting orbital 
entanglement in a mesoscopic system were proposed in normal-superconductmg[l] as well as purely normal 
systems [22, 23]. Energy-time entanglement [24] in a Franson like geometry [25] provides another spin- 
independent approach. In this paper, we discuss in detail the properties of the orbital entanglement and 
how it can be detected via violation of a Bell Inequality. In particular, we focus on the effect of decoherence 
on the entanglement. We follow closely the discussion in Ref. [1], but the results are of relevance also for the 
orbital entanglement studied in Refs. [22, 23]. An investigation of the effect of decoherence on entanglement, 
discussing the setup in Ref. [22], is presented by van Velsen et al in this issue. 
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2. Entangled two-particle wavefunction. 

We consider the system shown in Fig. 1. A single superconductor (S) [the upper and lower part in the 
figure connected e.g. via loop] is weakly coupled to a normal conductor, a ballistic two-dimensional electron 
gas, via two tunnel barriers 1 and 2 with equal transparencies r « 1. The normal conductor consists of 
four arms, 1A, IB, 2A and 2B, with equal lengths L. The arms 1A and 2A (IB and 2B) are crossed in 
a beam-splitter A(B) and then connected to normal reservoirs +A and —A (+B and — B). We consider 




Figure 1. The system of Ref. [1]: A single superconductor (S) is connected to four normal arms via two tunnel 
barriers 1 and 2 (thick black lines). The arms are joined pairwise in beam splitters A and B and end in normal 
reservoirs + and — . 

the low temperature limit, kT <C eV. A negative voltage — eV is applied to all the normal reservoirs and 
the superconductor is grounded. The voltage eV is smaller than the superconducting gap A, so no single 
particle transport takes place. It is further assumed that the size of the system is smaller than the phase 
breaking length. 

The superconductor as a source of spin-entangled electrons has been discussed in Refs. [11, 12, 16, 21]. 
As shown in Ref. [1], due to the coherent properties of the superconducting condensate, the superconductor 
can also act as an emitter of orbitally entangled pairs of electrons. To lowest order in tunnel barrier strength 
T <C 1, a pair of electrons is emitted either through contact 1 or contact 2, The state is thus a linear 
superposition of emitted pairs, an orbitally entangled two-particle state. In a second quantization language, 
the state of the entangled pair can be written as [1] 



l* out ) = |o) + |*) + |*), (l) 



where the states 



i*> = ^r^E e [^(^{(-^-^(^(-^jio), 

j=l,2r)=A,B 

l*> = T r dE Ha(E)c&{-E) c{ a (E)c\U-E) + 4 A (E)4U-E) 4 A (E)4U-E)} |0>. (2) 

4 J-eV 1 J 



Here, e.g. an operator c[ A (E) describes the creation of an electron with spin up at energy E, emitted 
through contact 1 and propagating towards beamsplitter A. The state |0) is the groundstate, a filled Fermi 
sea with electrons at all energies E < —eV and \^f) (|*)) describes orbitally entangled wavepacket-like states 
with two electrons propagating to different (the same) beamsplitters. As is shown in Ref. [1], the state |\&), 
where the two emitted electrons propagate towards the same beamsplitter, only contributes to the physical 
observables of interest to higher order in tunnel barrier transparency and can be neglected. In what follows 
below, we thus focus on the state \^f), describing one electron propagating towards beamsplitter A and one 
towards B. 
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To clearly visualize the orbital entanglement, we present the state |\E') in a first quantization notation. 



Introducing \1A, E}\ |) for the operator c\ A (E), the properly symmetrized wavefunction Y$>) is given by 



i*> 



reV 

/ dE {\IA, E)\1B, -E) + |LB, -E)\1A, E) + \2A, E)\2B, -E) + \2B, -E) \2A, E)) 

J-eV 



®(IT>U>-U>IT»- 



(3) 



Here the first and second ket correspond to the first and second particle respectively. We note that the 
total wavefunction is a direct product of an antisymmetric (the spin singlet state of the Cooperpair) spin 
wavefunction and a symmetric orbital wavefunctions. The orbital entanglement with respect to emission of 
the pair across barrier 1 or 2 can be even more explicitly seen by writing the wave function as 



I*) - |*12>®|*Afl>, |*12> = (|11> + |22))/A 



*ab) = f V dE(\A,E)\B,-E) + \B,-E)\A,E))®(\ T)| |) 

J-eV 



1)1 T)) 



(4) 



The state 1^12) is an orbitally entangled state with respect to emission across barrier 1 and 2. This two- 
dimensional space, called 12-space below, plays the role of a pseudo-spin space. The state I^ab) contains 
all additional information, such as energy and spin dependence. 

Propagating from the superconductors towards the beamsplitters the electrons pick up phases (fiA, </?ib, ¥>2A 
and ifi2B, where the index 1A etc. denotes the arm of the normal conductor, see Fig. 1. These phases can be 
taken energy independent for the system parameters considered. The beamsplitters are assumed to have an 
adjustable transparency and to support only one propagating mode. Such beamsplitters have been realized 
in recent experiments [26, 27]. The states | + ^4) and | — A) (| + B) and | — B)) for electrons going out 
into the normal reservoirs and the states \1A) and \2A) ( \1B) and |2£?)) of the electrons incident on the 
beasmsplitters are related via energy independent scattering matrices: 



A) 
A) 



COS <pA 

sin 4>a 



— sin q>A 

COS (f>A 



\2A) 
\1A) 



+ B) 
-B) 



COS (j)B 

sin 4>b 



- sin 4>b 

COS (j)B 



\2B) 
\1B) 



(5) 



Possible scattering phases of the beamsplitters are incorporated in the phases (fjrn with j = 1, 2 and rj = A, B. 
Under these conditions the beamsplitters only affect the state in 12-space, i.e. I^^)- On the normal reservoir 
side of the beamsplitters, the state ^12) can be written 



with the coefficients 

C+ + 



1*12) = c++| + +) + c__| - -) + c_+| - +) + c+_| + -> 



sin(^) s\n{4> B )e i{ipiA+vlB) + cos(0 B ) cos(4> A )e i(ip2A+V2B) /V2, 
cos(^)cos(0 i j)e l(¥,1A+VlB) +sm((f> B )sm((j>A)e l( - ip2A+V2B) ] /y/2, 



(6) 



cos(^a) sm(<j> B )e l{lplA+ipiB) - cos(0 B ) s\n{4> A )e i{v2A+V2B) 



c +- 



sm((p A ) cos(</> B )e 



sm{(j) B )cos((t)A)e i{ip2A+{fi2B) 



/V2, 
"2. 



(7) 



This provides a full description of the state of the particles exiting into the reservoirs. 



3. Detection via violation of a Bell Inequality. 

The entanglement gives rise to a nonlocal correlation between the two electrons, i.e. the outcome of a 
measurement at A (e.g. electron in A+) effects the probabilities for the outcome of the measurement in 
B (electron in B+ or B-). The implications of such non-local correlations have been at the heart of the 
discussion in quantum mechanics for almost a century [2, 3, 4]. Bell was the first to show [5] that the 
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predictions of quantum mechanics were incompatible with a whole class of theories, trying to explain the 
non-local correlations with local, realistic arguments in the spirit of Einstein, Podolsky and Rosen [2] . Bells 
proof was formulated in terms of an inequality which could not be violated by any system described by a 
local, realistic theory. Systems described by quantum mechanics could however violate such inequalities, 
showing that local realistic theories were untenable. 

Several approaches to Bell Inequalities (BI) in mesoscopic conductors have been discussed [28, 29, 30, 20] . 
Here we follow the approach in Ref. [1], which is closely related to Bells original derivation. Starting from 
the two particle state emitted from the superconductors, using the formulation in Eq. (6), we can directly 
calculate the probabilities for joint detection of an electron in reservoirs aA and one in j3B (a, (3 = +, — ), 
given by the simple relation 

P a p = \c a p\ 2 (8) 



This gives directly 



P++ = P— = [cos 2 (0a) cos 2 ((/) B ) + sin 2 (0 B ) sin 2 (0 j4 ) 

+ cos(2[ip A - (p B ])cos(<j) A ) sinO B ) cos(0 B ) sm(^ A )] /2 

P+- = P-+ = [cos 2 ((f>A)sm 2 ((f) B ) + cos 2 sin 2 (0^) 

- cos(2[<p,4 - ip B ]) co$,(4> A ) sinOs) cos(0 B ) sin(<^)] /2 



(9) 



where we have introduced the phases 2tfA — <Pia — ( P2A and 2tp B = tp 2 B — <Pib- The phases ipA and 
(fB can be varied independently and locally, at beamsplitters A and B, by e.g. electrostatic sidegates 
changing the length of the normal arms. Note that the joint detection probabilities are normalized such that 

P ++ + P-- + P-\ + P_ + = 1. Following Bell and later formulations [31, 32], we introduce the correlation 

function 

£(a,b) = P++(a,b) - P+_(a,b) - P_+(a,b) + P__(a,b). (10) 

where a (b) collectively denotes the settings of the phase <pA and angle 9a (phase <f>B and angle 9b)- In 
terms of these correlation functions for four different measurements settings, denoted a, a',b and b', the 
Bell Inequality can be formulated as 

-2<Sb<2, S B = £7(a,b)-£7(a,b') + S(a',b) + £7(a',b'). (11) 

where Sb is the Bell parameter. This provides us with a BI in terms of the joint detection probabilities P a p. 
However, in mesoscopic systems, in contrast to optics, these quantities are not experimentally accesible, they 
would demand a short-time coincidence measurement. The natural observables are instead zero-frequency 
current correlators, 



S a fi — 



/oo 
dt({5i aA {t)6i 0B (O) + Sif3B{0)6i aA {t)}} (12) 
-oo 



where SI ari (t) = I ari {t) — (I a ri) is the fluctuating part of the electrical current I ari {i) in reservoir aij. The 
zero-frequency current correlators are calculated along the lines of scattering theory [33] and are given by 
[1] 

4e 2 4e 2 
S ++ = S- = -^-\eV\P ++ = —\eV\P__ 

4e 2 4e 2 

S+-=S-+ = —\eV\ P+ _ = —\eV\P_ + (13) 

where P a p are just the joint detection probabilitcs of Eq. (9). Thus, the zero- frequency current correlators are 
directly proportional to the joint detection probabilitcs. The simple result can be understood by considering 
the properties of the time-dependent correlator (5I a A(t)Sl /jb(0)) . It decays as (r c /t) 2 for times t > r c , 
where t c = Ti/eV is the correlation time of the emitted pair. In the tunneling limit under consideration, 
r« 1, the correlation time is thus much smaller than the average time between the arrival of two pairs 
e/I <~ U/eVT 2 . As a result, only the two electrons within a pair are correlated with each other, while 
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electrons in different pairs are completely uncorrelated. Thus, the zero frequency current correlator in Eq. 
(12) is just a coincidence counting measurement running over a long time, collecting statistics over a large 
number of pairs. This leads to the important result that the Bell inequality, Eq. (11), can be directly 
formulated in terms of the zero-frequency current correlators in Eq. (12), i.e. the correlation function of Eq. 
(10) is given by 

E( a ,b) = [S++(a,b) - S+_(a,b) - S_+(a,b) + S__ (a,b)] /P (14) 

where P = (4e 2 /h)\eV\ is the proportionality constant between the current correlator and the joint detection 
probability. 

It is important to note that BFs in electrical conductors quite generally have to be investigated under 
conditions different from the ones in optics. In optics, the discussion has to large extent been focused 
on the detection problem, the fact that only a fraction of the photons emitted are actually detected [31, 
34, 35, 36] and how this effects the interpretation of an experimental violation of a BI. This is in general 
not a problem in electrical conductors, all electrons are detected in the reservoirs and contribute to the 
statistics. There are however different aspects of BI's to be considered in electrical conductors. First, 
the entagled quantities are not free electrons but quasiparticle excitations out of the Fermi sea. Second, 
there are several "loopholes" , due to possible interaction with phonons or other electrons or coupling to the 
electromagnetic environment. All these phenomena open up the possibility for local realistic theories that 
might explain the observed phenomena. In addition, more sophisticated local realistic theories allowing for 
some communication (possibly at the speed of light) between the detectors or the source and the detectors, 
the so called locality loophole, can hardly be ruled out in any mesoscopic experiment (but have in fact been 
done so experimentally [7] in optics). 

With these limitations in mind, a violation of a BI in an electrical conductor would neverthess provide 
a first indication of the nonlocal properties of correlations between electrons. More importantly, the BI 
treated completely within the framework of quantum mechanics can be used as a tool to extract the amount 
of entanglement. This becomes clear in the discussion below of the effect of dephasing on entanglement. To 
investigate the parameter dependence of the correlation functions in Eq. (14) we insert the expressions for 
the zero frequency correlators in Eq. (13) and obtain 

E(a, b) = cos(2(j) A ) coa(2(f> B ) + cos(2[^ - <p B ]) sin (20a) sin(20 s ). (15) 

We note that there are two types of parameters that govern these correlation functions, the phases ifA and 
tfB and the angles tf>A and 4>b ■ Putting the angles 4>a and 4>b to ir/4 (a semi-transparent beam-splitter) 
gives rise to the correlation function in terms of the phases ifA and <p b 

E(ipA, ipb) — cos(2[ipA - <Pb])- (16) 

Second, we can instead put the phases 2[ip A — <Pb] to a multiple of 2ir, giving a correlation function in terms 
of the beam-splitter angles 

E(<Pa,<I>b) = cos(20 A ) cos(20£>) + sin(20 A ) sin(20 B ) = cos(2[0 A - fo]). (17) 

This has exactly the same form as the correlation function in Eq. (16). Inserting these correlation functions 
into the Bell parameter in Eq. (11), and choosing e.g. <Pa = tt/8, 4>b — 7r/4, 4> a — 37r/8 and <p' B — n/2 (or 
equivalently the same values for the phases ipj V ) we get a Bell parameter S = 2\[2. This gives a maximal 
violation of the BI in Eq. (11). 

From this, one gets the impression that it is of no importance what parameters are adjusted. However, 
as investigated in detail below, the two approaches lead to very different results for the ability to violate the 
BI in the presence of dephasing. 

4. Effect of dephasing. 

There are many sources of dephasing in electrical conductors. Dephasing can arise due to direct in- 
teraction between different electrons or between electrons and phonons. Another source is coupling of the 
electrons to the electromagnetic environment. An investigation of dephasing in mesoscopic conductors is 
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given by Seelig, Pilgram and Biittiker in this issue. Focusing on systems with geometries and physical 
properties similar to the one in Fig. 1, we note that the effect of dephasing on the phase dependence of 
the conductance in Mach-Zehnder interferometers was recently studied in Refs. [37, 38]. Very recently, the 
effect of dephasing on the phase dependence of the conductance as well as the shot noise in Mach-Zcndcr 
interferometers was investigated experimentally [39] . Subsequently, a theory for the effect of dephasing on 
the shot noise in the same geometry was proposed [40]. We also note that the effect of dephasing on the 
possibility to violate a BI has been studied early in optics [41] . 

In this paper we study a simple model of the dephasing that nevertheless illustrates the general ef- 
fects of dephasing on the entanglement and the possibility to violate the BI. The dephasing is modelled 
phcnomenologically via a density matrix 

p = [|11)(11| + |22)(22| +7(|H)(22| + |H><22|)]/2 (18) 

which only affects the entanglement in 12-space. The dephasing parameter 7 can vary from 1, when the 
system is in a pure, coherent, fully entangled state and no dephasing is present, to 0, when the system is in 
a completely mixed state, fully dephased and with the entanglement suppressed. We note that this density 
matrix would e.g. arise from fluctuating phases ipj v , discussed in Refs. [37, 38, 40]. The density matrix in 
Eq. (18) gives rise to the correlation function, a modification of Eq. (15), 

E(a,b) =cos(2(/» A )cos(20 B )+7cos(2[^-<p B ])sin(20 A )sin(20 B ). (19) 

From this expression it is clear that the two approaches described above to adjust the parameters of the 
correlation function give very different results when trying to violate the BI. Fixing the angles <f)A and 4>b, 
the Eq. (19) becomes 

E(ip A ,ip B ) = 7COs(2[</m - tp B ]) (20) 

Thus, the dephasing parameter 7 takes on the role of what is known in the optical litterature as the visibility. 
For a visibility smaller than l/y/2, the BI can not be violated, which thus puts the rather strict condition 

7>^ (21) 

on the amount of dephasing allowed. Choosing the other approach above, fixing the phases ip A and ips, we 
find Eq. (19) modified to 

E((j> A ,(t> B ) = cos(2^ A )cos(2^ B )+7sin(2^ A )sin(2^ B ). (22) 

As was pointed out in Rcf. [1], this allows in principle for a violation of the BI for arbitrary dephasing. In 
the next section, this is investigated in detail. 



5. Optimal violation. 

We first try to find the angles which maximize the Bell parameter in Eq. (11). By maximizing the Bell 
parameter with respect to first (f> A and 4>' A and then to 4>b — <j>B> we arrive at the relations for the angles, 
all in the first quadrant 

tan(20 A ) = -7cot(0s), 
tan(2</>' A ) = 7tan(0s), 



where <fts — 4>b + §'b can be chosen at will. Inserting these angles into the Bell parameter, we arrive at 



S B = (24) 
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Figure 2. Left and middle figures, from Ref. [1]: The transmission probabilities Ta (dashed) ,T' A (dotted) and T' B 
(solid) as a function of Tb [T^ — cos 2 (0,,)], giving optimal violation of the Bell inequalities for dephasing parameters 
7 = 0.1 (left) and 7=1. (middle). Right figure: The fraction of the transmission probability space allowing for a 
violation, as a function of the dephasing parameter 7. 



showing that for an optimal choice of the angles <Pa,<I>'ai and 4>' B , the BI can be violated for arbitrarily 
strong dephasing, 7 > 0. The corresponding optimal transmission probabilities T v = cos 2 ((j> v ) are shown for 
7 = 1 and 7 = 0.1 in Fig. 2. The BI can thus in principle be violated for any amount of dephasing. However 
it might be difficult to produce beam splitters which can reach all transmission probabilities between and 
1. This is not a serious problem in the absence of dephasing, 7 = 1, a violation can be obtained for a large, 
order of unity, fraction of the "transmission probability space" (see Fig. 2). However, in the limit of strong 
dephasing, 7 <C 1, the set of transmission probabilities for optimal violation contains transmissions close to 
both and 1, see Fig. 2. Expecting unity transmission to be most complicated to reach experimentally, we 
note that by instead choosing transmission probabilities Ta = Tb = 0, T' B = 1/2 and T' A <C 7, the inequality 
in Eq. (11) becomes 2|1 + jT A \ < 2. This gives a violation, although not maximal, for all 7 <C 1. 

Apart from dephasing there are several other effects such as additional scattering phases, impurity scat- 
tering or asymmetric normal conductor arms that might influence the ability to violate the BI. These effects 
are discussed in Ref. [1]. Another possibility is that, due to asymmetries of the tunnel barriers Ti ^ T2, the 
amplitude for the process where the pair is emitted to |11) is different from the process where it is emitted 
to 1 22). This gives rise to a state, in 12-space, 

|*i 2 )= J_ (r 1 iii) + r 2 |22)). (25) 
v 1 1 + 1 2 

In this case [42, 43] we directly get the same expression for the coherence function in Eq. (19), with an 
effective 

21^2 

7 = TYTTT (26) 

This shows that it is in principle possible to violate BI for arbitrary asymmetry as well. 



6. Conclusions 

As a conclusion, we have investigated orbital entanglement in mesoscopic conductors with the focus on 
the effect of dephasing. The entanglement is detected via violation of a Bell Inequality formulated in terms 
of zero-frequency current correlations. We have investigated how the dephasing affects the possibility to 
violate the Bell Inequality and how system parameters can be adjusted for optimal violation. 
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